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Abstract. We report on a calculation of the growth of the mass of supermassive black holes at galactic centers 
from dark matter and Eddington - limited baryonic accretion. Assuming that dark matter halos are made of 
fermions and harbor compact degenerate Fermi balls of masses from 10 3 M© to 1O 6 M , we find that dark matter 
accretion can boost the mass of seed black holes from about ~ 5M© to 1O 3-4 M0 black holes, which then grow 
by Eddington - limited baryonic accretion to supermassive black holes of 1O 6 ~ 9 M0. We then show that the 
formation of the recently detected supermassive black hole of 3 x 1O 9 M0 at a redshift of z = 6.41 in the quasar 
SDSS Jl 14816.64+525150.3 could be understood if the black hole completely consumes the degenerate Fermi ball 
and then grows by Eddington - limited baryonic accretion. In the context of this model we constrain the dark 
matter particle masses to be within the range from 12 keV/c to about 450 keV/c . Finally we investigate the 
black hole growth dependence on the formation time of the seed BH and on the mass of the seed BH. We find 
that in order to fit the observed data point of Mbh ~ 3 x 1O 9 M0 and z ~ 6.41, dark matter accretion cannot 
start later than about 2 x 10 8 years and the seed BH cannot be greater than about 1O 4 M0. Our results are in full 
agreement with the WMAP observations that indicate that the first onset of star formation might have occurred 
at a redshift of z ~ 15 — 20. For other models of dark matter particle masses, corresponding constraints may be 
derived from the growth of black holes in the center of galaxies. 
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1. Introduction 

Over the past few years, the idea of dark matter (DM) and 
the possible existence of supermassive black holes (BH) of 
masses from 10 6 5 to 10 9 5 M Q at the center of galaxies 
(Macchetto et al. HW7\ Schodel et al. I2UU2|) have become 
deeply rooted within the astrophysical community. The 
exploration of the relationship between these two intrigu- 
ing problems of modern astrophysics has been the subject 
of an increasing number of papers. 

It has been established that the mass of the central 
BH is tightly correlated with the velocity dispersion a of 
its host bulge, where it is found that Mbh ~ cr 4-5 (Faber 
et al. \iWf\ Magorrian et al. IT§9"8l Ferrarese & Merritt 
I2TJUU1 Gebhardt et al. l2TMil Ferrarese l2TjTJ2l Haering & 
Rix 2004). This tight relation between the masses of the 
BHs and the gravitational potential well that hosts them 
suggests that the formation and evolution of supermassive 
BHs and the bulge of the parent galaxy may be closely re- 
lated, e.g. Wang, Biermann & Wandel l |2000|l . In addition, 
the recent discovery of high redshift quasars with z > 6 
(Fan et al. 2001 ) implies that the formation of supermas- 
sive BHs took place over fewer than 10 9 years. In spite of 
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the vast and tantalizing work undertaken on BHs, their 
genesis and evolution are not well understood (see Rees 
1984 for a review). Since the discovery of quasars in the 
early 1960s, it has been suggested that these objects are 
powered by accretion of gas onto the supermassive BHs of 
masses 10 6 — 10 9 M Q (Xvnden-Bell I196l?|) . Two scenarios 
have been discussed in modeling the growth of BHs. One 
is that BHs grow out of a low mass 'seed' BH through ac- 
cretion (Rees [BJ33JI; and another one is that BHs grow by 
merging (Barkana et al. 120011 Wang, Biermann & Wandel 
I2TJUU1 Gopal-Krishna, Biermann & Wiita EMl I2tM|> . In 
a recent paper, Duschl & Strittmatter (200|J have inves- 
tigated a model for the formation of supermassive black 
holes using a combination of merging and accretion mech- 
anisms. 

The purpose of this paper is to study the growth of 
BHs from dark matter and Eddington-limited baryonic ac- 
cretion. Also, we would like to obtain some constraints on 
the DM particle masses. In the past, self gravitating neu- 
trino matter has been suggested as a model for quasars, 

with neutrino masses in the range 0.2keV ~ to/ ~ 0.5MeV 
(Markov I1964fl . Later, neutrino matter was suggested 
to describe DM in clusters of galaxies and galactic ha- 
los with masses in the range of leV ~ to/ ~ 25eV 
(Cowsick & McClelland lUT751 Ruffini[T^). More recently, 
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fcrmion balls (FBs) made of degenerate fermionic matter 

of 10 kcV ~ Tfif ~ 25 keV were suggested as an alterna- 
tive to supermassive BHs in galaxies (Viollier 1994; Bilic, 
Munyaneza & Viollier 1993 Tsiklauri & Viollier 113351 
Munyaneza & Viollier I2002[l . It has been also suggested 
that if the Galaxy harbors a supermassive BH, then there 
should be a density spike in which dark matter (DM) 
falling towards the center could annihilate and the detec- 
tion of these annihilation signals could be used as a probe 
for the nature of DM (Bertone, Silk & Sigl l2UU21 Gondolo 
& Silk IT§9"9l and Merritt et al- HEHl . The current belief 
is that DM particles are bosonic and very massive, i.e. 

m dm ~ lGeV/c 2 . However, the absence of experimen- 
tal constraints on the weakly interacting massive particles 
(WIMPs) that probably constitute DM leaves the door 
open for further investigation of the hidden mass of the 
Universe. 

In this paper, we will assume DM to be of fermionic 
matter and described by a Fermi - Dirac distribution with 
an energy cutoff in phase space (King 196(J). We will then 
explore the limits for the DM particle masses in order to 
reproduce the mass distribution in the Galaxy (Wilkinson 
& Evans 11999(1 and then study the growth of a seed BH 
immersed at the center of the DM distribution in galaxies. 
We use degenerate FBs at the center of DM halos not as 
replacements for the BHs but as necessary ingredients to 
grow the BHs in galactic centers. 

The resulting distribution of stars around a massive 
BH was studied in detail in the 1970s and early 1980s 
in the context of globular clusters (Hills 119751 Frank 
fc Re es HWm Bahc all & Wo lf IT57H1 Duncan & Shapiro 
IT3S21 Shapiro EES). Peebles ijlHT^ studied the adiabatic 
growth of a BH in an isothermal sphere and showed that 
the BH would alter the matter density to an adiabatic 
cusp with p ~ r~ 3 / 2 . A few years later, Young (.1980) con- 
structed numerical models that confirmed Peebles' results 
and showed that the BH induces a tangential anisotropy 
in the velocity dispersion. In this paper, we use the results 
of previous calculations on the growth of BHs that are ac- 
creting stars (Lightman & Shapiro 1978 ) to investigate the 
growth of a BH that accretes DM. Here we assume that 
the physics driving the formation of the power law cusp in 
the star - star case is the same as in the case of DM par- 
ticle orbits being perturbed by molecular clouds. Julian 
(1967) investigated a similar scenario in which the stellar 
orbits in our Galaxy were perturbed by molecular clouds 
to explain the stellar velocity dispersion dependence on 
the star's age. Moreover, Duncan & Wheeler (1980) inves- 
tigated the anisotropy of the velocity dispersions of the 
stars around the BH in M87. 

Given the density distribution in DM halos, we are 
interested in establishing how a seed BH would grow by 
accreting DM. Moreover, the comparison of the growth 
of the BH from accretion of DM and Eddington - lim- 
ited baryonic matter would give us another piece of in- 
formation in the debate surrounding the nature of DM. 
We therefore investigate how a BH seed of 5 M Q could 



grow to a 3 x 10 9 M Q BH as recently detected in quasar 
SDSS J1148+5251 at z=6.41 (Willot, McLure & Jarvis 
2003J- Such a seed BH of a typical mass between 5 and 9 
Mq could have evolved in BH binaries (Podsiadlowski, 
Rappaport & Hau I2003[) . Here, we note that the seed 
BH could be an intermediate mass black hole (IMBH) 
of 10 3 ~ 4 M Q that might have formed from collisions in 
dense star-forming regions (Portegies Zwart & McMillan 
l2T)(T2l Coleman Miller IM)3l and van der Marel ISHnSI for 
a review). In fact, Wang & Biermann (1998) have estab- 
lished that a BH would grow exponentially with time ac- 
creting baryonic matter as long as the supply lasts. In 
addition, they were able to reproduce the observed cor- 
relation M BH /M sph using standard disk galaxy parame- 
ters. Assuming a cusp - like distribution of self-interacting 
DM (Spergel & Steinhardt lSUUUjl . Ostriker pUUUjl has es- 
timated that BHs could grow to 10 6 - 10 9 M Q from DM 
accretion. For completeness, we note that accretion of 
DM particles by BHs has recently been studied in Zhao, 
Haehnelt & Rees lj2TjU2)l and Read & Gilmore P7U5|> . 

Cosmological parameters of Hq = 70 km s _1 Mpc _1 , 
f^m = 0.3 and Q\ — 0.7 are assumed throughout this 
paper. In section 2, we establish the main equations to 
describe DM in galaxies. We then discuss the growth of 
seed BHs from DM and Eddington - limited baryonic mat- 
ter accretion in section 3 and conclude with a discussion 
in section 4. 



2. Dark matter in galaxies 

2.1. Main equations 

We characterize DM by its mass density p{r) and its veloc- 
ity dispersion a. DM is assumed to be collisionless and of 
fermionic matter and the mass of DM particles is denoted 
by mj. We will look for DM distributions with degener- 
ate cores at the center of the DM halos instead of a steep 
power law p ~ r~ 7 with 7 w 1 — 1.5 (e.g. Navarro, Frank & 
White 19 !J7f) as we are assuming here that a central cusp is 
a result of BH growth. In order to obtain bound astrophys- 
ical solutions at finite temperature, we follow King Q1966[) 
and introduce an energy cutoff in phase space using the 
Fermi-Dirac distribution function. We therefore adopt the 
following prescription for the distribution function 



f{E) = 8rf exp + 1 - - (1) 

[0 , E > E c , 

where gf is the spin degree of freedom, i.e. gf = 2 for 
Majorana and gf — 4 for Dirac's fermions, Ti is the Planck 
constant, E is the particle energy, p is the chemical po- 
tential and T stands for temperature. The constant C is 
chosen so that the distribution function vanishes at the 
cutoff energy E = E c . As a result, we get the follow- 
ing "lowered" Fermi-Dirac distribution that will be used 
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throughout the paper. 



f(E,E c ) 
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The mass density p(r) is given by 



p(r) 



4irf(p)p 2 dp(E) , 



(3) 



and the gravitational potential $ obeys Poisson's equation 

A$(r) = 4nGp(r) . (4) 

Next we introduce a dimensionless gravitational potential 
v, a radial variable x and an inverse temperature (3 as 



■($o-$(r)) 



a J 



(5) 



where a and To are the length and temperature scales, 
respectively. They are given as 



Zith A 



-i 2/3 



25/ 2 to4. 9/ G 3 /2M 1/2 



© J 



Tn = 



GMp 



Q m f 



ak 



(G) 



where k is Boltzmann's constant and the cutoff energy 
E c has been chosen as E c — m/$o- We also define the 
degeneracy parameter A as 



A = exp 



m/$ - p 
kT 



The equation for the gravitational potential can thus be 
written as 



Ifl = _^«-3/2_ 

xdx 2 2 P A 



A ( v(3 



with boundary conditions 



«(0) 



M 



BH 



Mr. 



v(x ) = , 



(8) 



(9) 



where xo — R/a is the dimensionless size of the DM halo 
at which the density and pressure of DM vanishes. It can 
be shown that for A — > oo, equation JSJ describes the 
isothermal gas sphere of King ( 1966) whereas for A — + oo, 
and v » 1 we recover the Lane-Emden equation for a 
degenerate FB fViollier Il994[) . J„(fc, rf) is a Fermi - like 
integral defined as 



rv i _ e £-77 

I n (A,r,)= / d£C- 
Jo 



1 + Ae^-i 



(10) 



The mass density p(r) can be written in terms of the di- 
mensionless potential v as follows: 



p(r) 



i)g f 2^ (GM { 



'0 \ a-3/2 A 



/3" 



A + l 



/2 



Using p(r) and the differential equation for the gravita- 
tional potential v, the mass enclosed within a radius r of 
the DM halo is given by 

E<E <(2) 

-r 

AiTp(r')r' 2 dr' 



M(r) 



= M Q (-xv'(x) + v(x) - u(0)) 



(12) 



and the total mass of the DM halo (including Mbh and 
the fermion halo) is given by 



M = M{R) = -Mqv'{x o )x , 



(13) 



where xq is the dimensionless size of the DM halo. We 
then study the solutions of equation JSJ for various values 
of A and mf. 

In Fig. 1, we plot the mass density and the mass en- 
closed for solutions with a 3 x 10 12 M Q mass and a radius 
of 200 kpc, which corresponds to our Galaxy. Throughout 
our paper, we will study solutions of type 3 which have 
degenerate FBs in the center. For the case of our galaxy, 
in order to have a degenerate FB of 3 x 10 6 Af Q at the 
center and fit the rotation curve of the Galaxy, i.e. repro- 
duce the 1/r 2 law, the mass of fermions should be about 
m/c 2 ~ 12keV. For a FB of 3 x 10 3 Af Q , one would need 
a fermion mass of ro/ ~ 450keV/c 2 . Thus under the as- 
sumption of a degenerate FB of 3 x 10 3 M© to 3 x 10 6 M Q 
at the center of a DM halo of 3 x 10 12 M Q with a density 
scaling as 1/r 2 in the outer edge of the halo, the fermion 



(7) mass is constrained in the range 



12keV~m f c 2 ~ 450 keV . 



(14) 



a,^].(h 

x 



Similar solutions for the DM distribution in our Galaxy 
have also been studied in Bilic et al. (2002a). 

The solutions shown in Fig. 1 have been obtained using 
the fact that at large distances the Fermi gas is non degen- 
erate. The degree of non-degeneracy is described by equa- 
tion (JJJ which contains the chemical potential p. Thus, 
the use of the chemical potential in the Fermi - Dirac 
ditribution allows to constrain the DM particles in the 
range from 12 to 450 keV, much below the Lee- Weinberg 
lower limit of about 2 GeV (Lee & Weinberg I1977|l which 
was established under the assumptions that DM parti- 
cles were in thermodynamic equilibrium at the freeze - 
out temperature. We also note that a suitable candidate 
for the DM particle as defined by the range given in l|14fl 
could be either the gravitino, postulated in supergravity 
theories with a mass in the ~ IkeV/c to - lOOGeV/c 
range (Lyth|2P0Q|, or the axino, with a mass in the range 
between ~ 10 and ~ lOOkeV/c , as predicted by super- 
symmetric extensions of the Peccei-Quinn solution to the 
strong CP problem (Goto & Yamaguchi I1992| ). For a re- 
cent review on DM particle candidates, the reader is re- 
ferred to Bertone, Hooper & Silk H2U04|) and Baltz 
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Fig. 1. The mass density and the mass enclosed within 
a radius r are plotted in the upper and lower panels, re- 
spectively. Near the center of the DM halo, the fermions 
are completely degenerate and condensed in a FB of mass 
- 3 x 10 6 A1©. In the outer edge of the DM halo, the den- 
sity p scales as 1 jr 2 (van Albada et al. I1985f) . The mass 
of fermions used in this plot is rrifc 2 w 12keV and the size 
of the FB is about 10 _2 pc. The data points for the mass 
within 50 kpc and 200 kpc are taken from Wilkinson & 
Evans l(T555j> . 

2.2. Black hole accretion rate and Paul! principle 

We define the local dynamical (orbital) time scale and the 
mass accretion flow as 



tdyn 5 
V ff 

Mbh = 4:TTr 2 p(r)vf f , 

where Vff is the free - fall velocity 

/2Gm(r)\ 1/2 
v ff = ~ 



(15) 
(16) 

(17) 



and m(r) is the mass enclosed within a radius r in the DM 
halo. 



Near the center of the DM halo, the core becomes com- 
pletely degenerate and the Pauli condition can be written 

as 

gj \ V 3 rrifV 

67?- 



" " ' 77,1/3 ^ 



(18) 



where n = p/mf is the number density for DM particles 
and v is the Fermi velocity. A seed BH of mass Mbh that 
settles at the center of the DM halo induces a free - fall 
velocity 



v ff 



2GM 



BH 



(19) 



Replacing v from the equation l|18|l with the last expres- 
sion for Vff, the density of DM at the Pauli limit is given 

by 

3/2 

// ■-- -r I I • (20) 



67t 2 h 3 \ r 

We then rewrite equation (|16f) for the mass flow using the 
Pauli flow condition near the BH 



M BH = 



m)G 2 M 2 



8g/ 

1.03 x 10~ 7 



BH 



fJf 



( M 



BH 



rrifc 
15 keV ) \ Ma 



M yrt21) 



Thus as soon as the region of BH influence is reached, 
the accretion flow becomes independent of the radius. 
However, the BH accretion rate depends on the fourth 
power of the fermion mass and on the square of the BH 
mass. This dependence is of course a consequence of the 
Pauli principle applied near the BH. 

Here we note that the mass Mf of nonrelativistic de- 
generate FBs scales with radius Rf as 



M F Rt 



/15keV\ 
27.836M© ( — — 



V m f c2 ) 



— PC 



(22) 



The physics of degenerate FBs as a model of DM at the 
center of galaxies was considered in a series of papers (Bilic 
et al. I2002al20()2bl Bilic, Munyaneza & Viollier 
Munyaneza, Tsiklauri & Viollier lT^gllTMal Munyaneza & 
ViollierEOni Tsiklauri & Viollier ITMBI and Viollier IT§9l|) . 
In these papers, it was argued that the supermassive com- 
pact objects in galactic centers are FBs rather than BHs. 
Our approach in this paper is that FBs are not a replace- 
ment for BHs but are a necessary step for them to grow 
in galactic centers. 

In Fig. 2 we plot the mass - radius relation for FBs 
for different values of the fermion mass. In the same plot 
we have drawn the BH line and two horizontal lines that 
correspond to a mass of 3 x 10 3 M Q and 3 x 10 6 Af Q . From 
this plot, we find an upper limit for the fermion mass of 



to/c 2 < 32 MeV 



(23) 



This upper limit is the value of the fermion mass that 
would correspond to the point in the graph where the BH 
line crosses the horizontal line of a mass Mf = 3 x 1O 3 M0. 
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Fig. 2. The mass - radius relation for degenerate FBs. The 
total mass scales as Mp ~ Rp 3 mJ 8 . In the same plot, we 
have shown the BH line and two horizontal lines for the 
lower and upper limits for the total mass of the degener- 
ate Fermi balls. Relativistic effects make the mass - radius 
relation curve have a maximum mass Mov (see equation 
124(1 and the size of such a FB at the Oppenheimer Volkoff 
limit is only AAhR s ov where Rqv ^ s the Schwarzchild ra- 
dius of the mass Mov- The line of BHs and a FB at the 
OV limit are nearly indistinguishable in this plot. 



The choice of the lowest value for the FB mass is mo- 
tivated by the fact that once the BH has consumed a 
FB of 3 x 10 3 M Q , it would grow by Eddington - limited 
baryonic matter accretion to 3 x 10 9 M Q as observed in 
quasars at redshifts z ~ 6. The choice of the second value 
of 3 x 10 6 M Q , is due to increasing evidence for the ex- 
istence of a BH of mass 3 x 10 6 Af Q at the center of our 
galaxy. Thus we will assume throughout this paper that 
galaxies which have BHs had degenerate FBs of masses 
ranging from ~ 10 3 A/q to ~ 10 6 Af Q . The fermion mass 
range given by equation ((14(1 is of course well below the 
upper limit for the fermion mass from equation 1(23(1 . 

For completeness we note here that stable degener- 
ate FBs exist up to a maximum mass also called the 
Oppenheimer Volkoff limit 



A I, 



ov 



0M195M Pl mJ 2 g f 1/2 

= 2.7821 x 10 9 A/ Q ( i^XV (24) 
\ m fC J \9fJ 

where M Pt = (hc/G) 1/2 is the Planck mass (Bilic, 
Munyaneza & Viollier 11999(1 . The size of a such FB at 
the Oppenheimer Volkoff limit is only Rov = 4.45i?Qy 
where R ov is the Schwarzschild radius of the mass Mov- 
Thus there is practically no difference between a black 
hole and a FB of the same mass at the OV limit, as the 
closest stable orbit around a non-rotating BH has a radius 
of 3R S only. 

From equation 124(1 . the fermion mass needed to form 
a degenerate FB of 3 x 10 9 Mq at the OV limit is about 
rrif ~ 12keV/c 2 . This value is in agreement with the lower 
limit on the mass of fermions obtained from fitting the 



Fig. 3. Fit of the rotation curve of dwarf galaxies DD0154. 
The data points can be fitted by a non-degenerate Fermi 
- Dirac distribution function for DM particles in the mass 

range discussed in this paper, i.e. for 12keV ~ rrifc 2 ~ 
450keV. In this plot we have used a fermion mass of m,j ~ 
305/ceV/c 2 to fit the rotation curve of the dwarf galaxy 
DD0154. Here we also note that some small addition of 
baryonic matter would probably allow a closer fit to the 
data. 

Galactic DM halo with a FB of 3 x 1O 6 M at the center 
and 1/r 2 density law from about 3 kpc. The most compact 
supermassive object of 3 x 1O 9 M0 has been identified as 
a BH at the center of M87 (Macchetto et al. 11997(1 . 

Apart from fermions as candidates of DM, bosons have 
also been considered as an alternative to the DM con- 
stituents in galaxies. Self - gravitating bosons can form 
boson stars and they are prevented from complete gravi- 
tational collapse by the Heisenberg uncertainty principle 
ArAp ~ Ti (Schunck & Liddle HU77|l . Boson stars are de- 
scribed by Einstein's field equations coupled with the non- 
linear Klein Gordon equation for a complex field ^ with 
self - interaction of the form U(\^>\) = m 2 B \^>\ 2 + re|*| 4 
where is the boson mass and k is a coupling constant. 
It has been shown that boson stars have a maximum mass 
called the Kaup limit given by 



M. 



Kaup 



0.633M Pl /m B 



(25) 



for k = 0. Thus for a boson mass of ~ 10 _23 GeV/c 2 
the Kaup mass limit is about 1O 6 M (Torres, Capozzicllo 
& Lambiase 120001 . For ms~ lGeV/c 2 the Kaup limit is 
as small as 10 _17 Mq. However if k ^ 0, the combination 
of these two free parameters ms and n let boson stars have 
any mass for any chosen value of the boson mass tub - It is 
in this framework that boson stars were attractive as they 
could form very massive objects with mg > lGeV/c 2 ( 
Schunck & Mieke I2UU5}> . 

In Fig. 3, we fit the rotation curve of dwarf galax- 
ies using a Fermi - Dirac distribution. The rotation curve 
data points are well fitted in the outer region. In the in- 
ner region, i.e. around 1 kpc, baryonic matter is perhaps 
needed to get a better fit of the rotation curve. In order to 
fit the rotation curve of dwarf galaxies, a non-degenerate 



6 



F. Munyaneza & P.L. Biermann: Growth of black holes 



FD distribution is needed. The data points are taken from 
Carignan & Purton (fTggBf and the mass of the the dwarf 
galaxy DD0154 is taken to be 3 x 10 9 M Q with a size of 
~ 10 kpc (Kravtsov et al. I1998[) . The merging of dwarf 
galaxies would generate another galaxy with a degenerate 
FB, and the BH would grow by the same mechanism con- 
sidered in this paper in the merged galaxy. The Fermi - 
Dirac distribution studied in this paper allows two types 
of solutions: One which fits dwarf galaxies with a non - de- 
generate core and the other solution with 1 /r 2 drop off in 
the outer edge of the halo. Degenerate FBs are consistent 
with data on the rotational curves only in the second type 
of solution with a 1/r 2 density profile at large distances. 

3. Black hole growth 

3.1. Inner dark matter accretion and the quantum 
cascade mechanism 

In this section we will consider the accretion of inner DM, 
i.e. DM inside the FB and in the next section we will anal- 
yse the accretion of DM outside the FB. In the last section, 
we established that the inner part of the the Fermi ball is 
completely degenerate and have used the Pauli Principle 
to obtain the BH consumption rate given by equation H21|) , 
which can be integrated to give the BH mass as a function 
of time: 



M BH = M Q 



Mo 



M 



BH 



1 



t-e 



(26) 



where M l BH is the mass of the seed BH at the starting 
time t* of accretion, and to is a time parameter given by 



t = 9.7 x 10V ( ^) g f 



(27) 



The last two equations define the growth of a seed BH 
immersed at the center of a DM halo with a degenerate FB 
from which it is seen that the BH mass strongly depends 
on the fermion mass. 

From equation (|26|) . we can estimate the time too 
needed for the BH mass to be infinite to be 



to 



M BH 



-to 



(28) 



which for an initial seed of 5M Q at t % = would give 



1.94 x 10V 



/15keVV 



\ mfc 2 



(29) 



For a fermion mass of rrif = IkeV/c the above timescale 
too would be about 5 x 10 10 years whereas a fermion mass 
of rrif = lMeV/c 2 would give a growth timescale of about 
too = 5 x 10~ 2 years. However, the BH mass cannot be 
greater than the mass of the FB that feeds it and we will 
require that the BH mass will grow from DM accretion 
only up to the mass of the FB and this will constitute 
one of the main assumptions in this paper. Once the BH 
mass has grown to the FB ball mass, a further growth 



to 10 6 9 Mq is achieved via Eddington - limited baryonic 
matter accretion. 

The mechanism of consumption of DM particles by the 
BH can be understood using a quantum cascade mecha- 
nism. First, low angular momentum DM particles at the 
inner Fermi surface in phase space will be consumed by 
the BH. Due to a high degeneracy pressure, high angu- 
lar momentum particles will be pushed to the inner Fermi 
surface and taken up by the BH. Outside the FB, fermions 
will hit the outer Fermi surface and for the FB to be hit 
continuously by DM particles, we use molecular clouds as 
perturbers of DM orbits. Thus, we have a three body in- 
teraction between the FB, the molecular clouds and DM 
particles. For our mechanism of consumption of DM by 
the BH to work, we will assume that the dark matter par- 
ticles interact with the FB through inelastic collisions. In 
hitting the FB, some of the DM particles will be excited to 
a slightly higher energy level and get stuck in high momen- 
tum segments of phase space. This process will continue 
until the FB has been completely consumed by the BH 
and thereafter baryonic matter accretion will control the 
BH growth. 



3.2. Outer dark matter accretion 

We study the growth of a seed BH of about 5M Q in- 
troduced at the center of DM halos of masses ranging 
from 1O 9 M to 10 14 ' 5 M Q . The lower limit corresponds 
to dwarf galaxy halo masses while the upper limit corre- 
sponds to the most massive halo such as that of M87. The 
gravitational potential and the mass density of DM are 
given by Eqs. J8J and JTIJ, respectively. It is well estab- 
lished that the BH will induce a density law p(r) ~ r~ 3 / 2 
at the center of a DM halo (Lightman & Shapiro 119781 
Quinlan et al. lT535l MacMillan & HenriksenEIEl Merritt 
I2UU3|| . Following Spitzer & Schwarzschild ^1951111953( 1 and 
Julian l|1967[l . we assume that DM orbits outside the de- 
generate FB are perturbed by molecular clouds of mass 
mpert ~ 10 6 M©. It has been recently reported that molec- 
ular gas of total mass ~ 10 10 Af© is observed in the host 
galaxy of the same quasar SDSS J114816.64+525150.3 at 
z w 6.41 (Walter et al.E003: Carilli et al. 



In this paper, we assume that there is a compact de- 
generate FB at the center of galaxies that have BHs. The 
seed BH first consumes the degenerate FB and grows to 
a mass of about 10 3 ~ 6 Mq and then grows by Eddington 
- limited baryonic matter accretion to higher masses of 
10 6 - 9 M Q . 

We will be mostly concerned with the calculation of 
the time needed to refill the BH loss cone using molecular 
clouds as DM perturbers. To this effect, we denote by r the 
radial distance from the BH to the DM particle position, 
f the distance from the BH to the molecular clouds, 9 is 
the angle between the direction from the BH to the DM 
particle and its projection in the disk and ip is the angle 
between the direction from the BH to the molecular clouds 
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Fig. 4. Geometry for the system BH (BH), DM (DM) par- 
ticles and molecular clouds (MC). FD stands for the foot 
(projection) of the DM particle in the disk. The angles BH 
- FD - DM and MC - FD - DM are both right 
angles. 



and to the foot of the DM particle location in the disk. 
The diagram for the geometry used is shown in Fig. 4. 

In order to get the expression for the time T re fui(r) 
needed to refill the loss cone which is due to multiple small- 
angle gravitational (Coulomb) encounters, we evaluate the 
random walk integral (Spitzer 1987) 



Av = 



( 2vpo 
V P 



Fig. 5. The refilling and orbital timescales plotted as a 
function of the radius from the center of the DM halo. 
In this plot, the mass of the FB is 3 x 1O 6 M and the 
corresponding fermion mass is 12keV/c 2 . For small r, the 
time to refill the loss cone is much greater than the age 
of the Universe. For the growth of the BH, we use the 
orbital time scale and the BH feeds on DM from inside 
the FB. For the DM inside the FB to be fully consumed 
by the BH, wc introduce a quantum cascade effect which 
is a result of the existence of a high degeneracy pressure 
that pushes high angular momentum particles to the inner 
Fermi surface. The particles at the inner Fermi surface in 
6D phase space will then be consumed by the BH. This 
process continues until the BH has consumed the entire 
FB. The BH then grows by Eddington -limited baryonic 
matter accretion. The FB is fed from DM from outer orbits 
but its accretion rate is much lower than that of the BH 
inside the FB as seen in the next figure. The vertical thick 
line shows the size of the FB. 

density into equation H30|) . after integration we get the 
following equation for the time to refill the loss cone: 



rdrdipn pert (r)vT re fui(r) ,(30) 



efui(r) 



v 3 r 



where v is the relative velocity of DM particle and molec- 
ular clouds at large separation, po and p are given by 



2ttG 2 (m DM + m pert ) a I(r,9) \ v 



Av 



with 



Po = 



G(m pert + m DM ) 



P 2 = 



-r 2 



—2rf cos 9 cos ip 



fAv 
•(31)(V 



R, 



( 2GM BH 
\ c 2 r 



(33) 



(34) 



In the above equations, we have used the values of 
0.1 kpc and 20 kpc for r m in and r max , respectively. We 
will also assume the following 3D number density distribu- 
tion equivalent to the surface density for molecular clouds 
(Glisten & Mezger 11352)1 



where the function I(r, 9) is given by the following integral 

exp(— x) 



I(r, 9) 



dx- 



2(£)acos(20)z3 + (£)4 



> (35) 



n P ert(r) = — exp 

2r V r 



(32) 



where rg = 5 kpc is a scaling distance and is a con- 
stant denoting the surface density to be chosen a pos- 
teriori. Inserting the last expression (|32|l for the number 



In the last integral, x min = r mm /r = 0.02 and x max = 
?max/ r o = 4. The function I(r,6) tends to a constant for 
7-^0 whereas for large values of r it scales as 1/r 2 . 

The FB is consumed from inside by the BH with an 
accretion rate 



Mbh = 4:irr 2 pvff, r < Rp, 



(36) 
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Fig. 6. The mass accretion flow as a function of distance 
from the center. For a given mass of the BH, the accretion 
rate is independent of the radius until the Bondi radius 
where the BH dominates the gravitational potential of the 
fcrmions. As the BH mass approaches that of the Fermion 
ball, the BH accretion rate becomes constant until the 
radius of the FB and as the size of the FB is reached, the 
DM accretion rate onto the FB drops dramatically. Thus 
the growth rate of the FB is much lower than that of the 
BH which grows by feeding on DM from inside the FB. 
Similarly to the last plot, we have used here a FB mass of 
3 x 10 6 M Q . 

The left side of the plot shows the BH accretion rate and the 
FB consumption rate is shown on the right side of the figure. 

which after using the Pauli degeneracy condition is given 
by equation l|21[) . Thus the BH gains by consuming DM 
inside the FB and of course the FB loses its DM parti- 
cles. On the other hand the FB is also fed on the outside 
by direct collision via the loss cone mechanism. The FB 
growth rate is given by 



M 



FB 



Anr 2 p{r)8VL{r)T{r)~ l dr r > R F , (37) 



where Rp is the size of the FB, Rdm is the DM halo size 
and Sn, is the fraction of the solid angle subtended by the 
loss cone and is given by 



5CI = 7T — 



( 2GMbh 
\ c 2 r 



The timescale r for the FB growth is given by 

^dyn ) ^refill ^ ^dyn 



(38) 



(39) 



refill 



; tdyn ^ ^refill ; 



where td yn and t re fm are the dynamical (orbital) and loss 
cone refilling timescales, respectively. As soon as the FB 
is consumed, the capture radius for DM particles coming 
from the outside drops rapidly. 

In Fig. 5, we plot the dynamical and the refilling 
timescales as functions of the radius. At small radii, i.e. 



Fig. 7. The black hole mass growth from only DM accre- 
tion. The fermion mass is varied as shown on the plot. 
It can be seen that the BH mass can grow to 10 3_4 Mq 
in about 10 7-8 years if the fermion mass is in the range 
between about 1 keV and 1 MeV. For fermion masses 
rrif ~ IkeV/c 2 , the formation of a BH of 10 3 M© would 
require a timescale much greater than than a Hubble time. 

On the other hand, for fermion masses m/ ~ lMeV/c , 
a 10 3 M© BH would have formed in less than than one 
year after the beginning of the Universe. Thus the de- 
generacy pressure of fermions allows for the capture of a 
large amount of mass if the fermion mass is constrained 
in a range between about 1 keV and 1 MeV. The data 
point shown at t = 8.47 x 10 8 years is the recent detec- 
tion of a BH mass of 3 x 10 9 M Q BH in the quasar SDSS 
SDSSJ114816.64+525150.3. 



r < R s , the time t re fui to refill the loss cone is too long for 
the DM to feed the BH. Near the BH, we use the dynam- 
ical timescale to feed the BH. For radii greater than the 
size of the FB, we will use the slowest of the two timescales 
to feed the FB. Thus the FB feeds from the outside DM 
particles via the loss cone mechanism while the BH feeds 
from DM inside the FB via the quantum cascade mecha- 
nism. 

In Fig. 6, we plot the mass flow as a function of the ra- 
dius from the center. To this end, we fix the mass of the FB 
to 3 x 1O 6 M and use a fermion mass of rrif ~ 12keV/c . 
Using the BH accretion rate equation H21|) . the accretion 
rate becomes M ~ lO 6 M0yr -1 for a fermion mass of 
rrif = 12keV/c 2 andaBH mass of 3x 10 6 M Q . For fermions 
of mass 450 keV/c 2 and a BH mass of 3 x 10 3 M Q , the ac- 
cretion rate is M ~ 10 4 M Q yr _1 . It is seen than the BH 
accretion rate is independent on the radius and drops sig- 
nificantly when the size of the FB is reached. Thus the BH 
grows much faster from the inner DM particles than the 
FB from the outer DM particles. After the BH has reached 
a mass equal to that of the FB, i.e. 1O 3_4 M0, its growth 
is controlled by Eddington baryonic matter accretion. 

In Fig. 7, we plot the BH mass as a function of time 
for different values of the fermion mass. The BH mass - 
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Fig. 8. The mass Mbh of the BH is plotted as a function 
of time t for DM and Eddington - limited baryonic mat- 
ter accretion. It is seen from this plot that DM dominates 
the BH growth at early times. The mass of fermions is 
varied as shown on the plot. The total mass of the degen- 
erate FB has a mass of 3 x 1O 3 M0. We find that fermion 

masses of my ~ 4keV/c suffice for stellar seed BHs to 
grow to about 1O 3 M which then grow by Eddington - 
limited baryonic matter accretion to about 3 x 1O 9 M . 
The obtained constraint on the fermion mass covers the 
range of dark matter masses discussed earlier. The data 
point shown at z = 6.41 corresponds to the BH mass of 
- 3 x 1O 9 M in the quasar SDSSJ114816.64+525150.3. 



time dependence is given by equation (|26l) . It is shown 
that for fermion masses lower than about 1 keV, it would 
take more than 10 10 years for the BH to grow to 1O 3_4 M . 
However, due to the Pauli principle used in the derivation 
of equation l(26|l , if the fermion mass is constrained in the 
mass range between about 1 keV and 1 MeV, a seed BH 
of 5 Mq can grow to a 1O 3_4 M BH in about 10 7 " 8 years, 
which would then grow by Eddington - limited baryonic 
accretion to lO 9 Af at a redshift of z w 6.41 (see next 
sub - section) . Fermion masses greater than about 1 MeV 
would make the BH grow much faster and the BH could 
be as large as 1O 3_4 M in less than about one year after 
the beginning of the Universe. 

In the next section, we will discuss the Eddington bary- 
onic matter accretion and discuss the plots of BH growth 
from both DM and Eddington baryonic matter accretion. 

3.3. Eddington limited baryonic accretion 

The details of Eddington - limited baryonic accretion have 
been studied by many authors, here we refer the reader 
to Wang & Biermann (1998) and provide just the main 
formulae: 

dMsH = M BH ^ 
dt tedd 



Fig. 9. The same as in the last figure but for a degenerate 
FB mass of 3 x 1O 6 M . 

where tedd = 10 769 years x (jfj) is the corresponding 
time scale, e being the efficiency. This value for tedd has 
been chosen to fit the exponential growth of BHs given by 
Fig. 1 of Wang & Biermann 1)1998(1 . The BH mass is found 
to grow in this case as 

M BH (t) = M BH exp (^-) , (41) 
\teddJ 

and the corresponding Eddington luminosity for an effi- 
ciency e of 0.1 is 

L , e M M c 2 =0.1^c 2 . (42) 
tedd 

In Fig. 8 and 9, we present the growth of supermassive 
BHs accreting DM and baryonic matter. It is seen from 
the plot that a BH grows faster in the beginning accreting 
DM rather than baryonic matter. In Fig. 8, the degenerate 
FB has a mass of 3 x lO 3 Af whereas we use a FB mass 
of 3 x lO 6 Af in Fig. 9. The data point shown in Fig. 8 
& 9 at t = 8.4 x 10 8 years is the recent detection of a 
supermassive BH of mass of Mbh ~ 3 x 1O 9 M in the 
quasar SDSS J114816.64+5251 at a redshift of z = 6.41 
(Willott, McLure & Jarvis l2TiTl3l Fan et al. EUUT)) . 

It is seen from Fig. 8 & 9 that only baryonic matter 
accretion cannot fit this data point and the curve for the 
growth of the BH due to baryonic matter lies on the right 
side of the data point. Using equation (|41|l the total mass 
of the BH from Eddington - limited baryonic accretion at 
t = 8.4 x 10 8 years i.e. z = 6.41 is 1.4 x 10 8 M o which 
is about ten times less than the inferred BH mass of 3 x 
lO 9 Af , clearly outside the error range. The growth of a 
seed BH from 5 Mq to 3 x 1O 9 M is accomplished in two 
steps: First, the BH completely consumes the degenerate 
Fermi core with masses between 3 x 1O 3 M and 3 x 1O 6 M . 
In the second step, the 10 3 — 1O 6 M BH feeds on baryonic 
matter accretion to reach the mass of 3 x 1O 9 M at a 
redshift of z ~ 6.41. 

In Fig. 10, we investigate how the growth of the mass 
of the BH depends on the accretion starting time. To this 
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Fig. 10. The BH growth dependence on the starting time 
of accretion. We find that starting DM accretion at a time 
of 2 x 10 8 years would fit the recent discovery of a 3 x 10 9 Mq 
BH in the quasar SDSSJ1148+5251 at a redshift of z = 
6.41, which is the most distant known quasar, observed 
only 840 million years after the beginning of the Universe 
(Willot, McLure & .Iarvis l2T)()3[) . The fermion mass used in 
this plot is rrif = lOkeV/c . However, any higher fermion 
mass would more sharply and rapidly raise the mass of a 
seed BH from 5 M to ~ 10 3 M©. 



effect, we fix the total mass of the FB to 3 x 10 3 M Q and the 
fermion mass to m/ = 10keV/c 2 and vary the initial time 
when the DM accretion starts as shown on the graph. We 
find that the data point of 3 x 10 9 M Q at z ~ 6.41 is well 
fitted for an accretion starting time of t ~ 2 x 10 8 years 

> 2 

for any fermion greater than m/ ~ lOkeV/c . Here we 
also note that Eddington baryonic matter accretion could 
also allow the seed BH to grow to a 3 x 10 9 Af Q BH if the 
efficiency is e = 1.1 x 10~ 2 . The accretion starting time of 
N2x 10 8 years corresponds to a redshift of reionization 
of z w 17 as obtained by WMAP observations (Krauss 

Finally in Fig. 11, we investigate how the growth of 
the BH is affected by the mass of the seed BH at different 
starting times of accretion. In order to constrain the mass 
of the seed BH, we vary the mass of the FB from 10 3 A/ Q 
to 10 6 M Q and also the mass of the seed BH from 5M Q 
to 1O 4 M at initial times of t l = 3.5 x 10 6 and 2 x 10 s 
years. These values correspond to the earliest and latest 
plausible times for the formation of the first stellar seed 
BHs in order to grow into supermassive BHs at redshift 
z w 6.41. From recent WMAP observations (Spergel et 
al. I2003[l . BHs of masses ~ 10 5 Af Q might have formed 
at a redshift z > 20. In Fig. 11, we have plotted only the 
limiting curves which show that the mass of the FB should 
be Mp < 1O 4 M and the mass of the seed BH cannot be 
greater than about 1O 4 M . The best fit is obtained for a 
BH formation time a reionization time oft 1 ~ 2 x 10 8 years, 
i.e. z ~ 17. 



Fig. 11. We investigate how the BH growth depends on 
the mass of the seed BH. To fit the data point of mass M ~ 
3 x 10 9 M Q at z w 6.41, the seed BH cannot be greater than 
10 4 Mq for a FB of 10 4 M Q and the accretion starting time 
should be about z m 17. In this plot, the fermion mass is 
fixed torn/ ~ 10keV/c 2 . Thus, any other mechanism that 
could create massive seed BHs would generate the same 
growth. 



It is interesting to estimate the maximum luminosity 
attainable by an accreting BH of this mass. This is usually 
considered to be the Eddington luminosity at which the 
outward radiation pressure equals the inward gravitational 
attraction. Using equation H42fl . the luminosity of a 3 x 
10 9 A/ Q BH is of the order of a few 10 47 erg s" 1 and the 
corresponding accretion rate is about 65 Af Q /yr. However, 
such high luminosity values (Fan et al. I2UU1|) would be 
obtained only for a standard efficiency of e = 0.1 which as 
seen from Fig. 8 & 9 does not fit the data point for only 
baryonic matter accretion. 

3.4. BH masses from energy conservation arguments 

A gravitational system in which molecular clouds provide 
the interaction path between DM particles and feeding the 
BH leads to a strong limit on the maximum mass of the 
BH. For an initial distribution of DM in the halo without 
a BH, the total energy of the system is given by 



Ei = 



5R D 



M 



(43) 



where M^m, Rdm stand for the total mass of the DM 
halo and its size, respectively. After the BH is formed at 
the DM halo center, its total energy is given by 



E 



f 



-ZGM% H 
5i? s 



with R s being the Schwarzschild radius 

_ 2GMbh 
n s = . 



(44) 



(45) 
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Equating Ef to Ei as a limit, we arrive at the following 
constraint 



M. 



DM 



f Rdm\ 



1/2 



M BH \ Rs J 

from which wc derive the BH mass limit 



Mmax 
BH 



c 2 Rdm 



(46) 



(47) 



Assuming an isothermal gas sphere distribution for the 
DM halo, i.e. M DM = 2cj 2 Rdm/G, with a being the DM 
velocity dispersion, we arrive at the final expression for 
the BH mass 



Mmax 
BH 



or 



Scr Rdm 

c 2 G 



M max = 4 ( ~_ 



M, 



DM 



(48) 



(49) 



The last two equations constrain the properties of the DM 
matter halo to the BH and could be considered as the the- 
oretical interpretation of the Mbh — cr observed correla- 
tion in elliptical galaxies. Here, we note that a stands for 
the velocity dispersion of DM. Using the above equations, 
the BH mass limit Mgg* is obtained to be ~ 10 2 Af Q , 
~ 10 6 M o , and ~ 10 9 Af Q for dwarf galaxies, the Milky 
Way and M87, respectively. The corresponding DM ha- 
los have masses of ~ 4 x 10 9 Af Q , ~ 3 x 1O 12 M and 
~ 4 x 10 14 M Q . In order to get the above estimates for 
the maximum BH mass, we have used a DM halo size of 
Rdm ~ 9 kpc for dwarf galaxies, Rdm ~ 200 kpc for the 
Milky Way and the size of M87 is assumed to be ~ 5 Mpc. 



4. Discussion & conclusion 

In this paper, we have investigated the growth of a stel- 
lar seed BH immersed at the center of DM halos with 
degenerate FBs of mass from ~ 10 3 Af Q to ~ 10 6 Af Q . 
Using the Pauli exclusion principle, we have established 
that the BH accretion rate strongly depends on the mass 
TO/ of the fermions as Mbh ~ mjAf^jj and thus estab- 
lish for the first time the relationship between BH growth 
and fermionic DM. We have shown that in order to fit 
the DM distribution in the Galaxy with such degener- 

2 <C 

ate cores, the DM particles should be in the 12keV/c ~ 

TO/ ~ 450keV/c mass range. We have shown that such 
DM masses could be used to fit the distribution of DM in 
dwarf galaxies. FBs of masses 10 3 ~ 6 Af could only ex- 
ist in such galaxies where the density drops off as 1/ r 2 at 
large distances. Dwarf galaxies as well as cluster of galax- 
ies do not host FBs as their data can only be fitted by a 
non - degenerate Fermi - Dirac distribution of the King 
type. We have argued that the merging of dwarf galaxies 
would lead to the formation of galaxies with degenerate 
FBs. The wideness of the fermion mass in the keV range 
is due to the FB mass range from 10 3 to 10 6 M Q that we 
have adopted in this paper. Our main assumption is the 



use of fermions as DM candidates in galaxies. However, if 
one uses bosons instead of fermions, the range of DM par- 
ticle masses would of course differ from the one obtained 
in this paper. The range of fermions used in our paper is in 
conflict with the Lee- Weinberg lower limit of ~ 2GeV/c 2 
on the fermion mass (Lee & Weinberg 1977). This is due to 
the fact that the derivation of the Lee- Weinberg limit as- 
sumes a freezout from equilibrium distributions. Our sim- 
ple model uses the chemical potential which allows for non 
equilibrium distributions and this might modify the Lee- 
Weinberg argument; this remains to be demonstrated. 

If we use heavy fermions with masses of about 
lGeV/c 2 , then according to equation (|24H the maximum 
mass allowed for the degenerate Fermi core would only 
be of about lAf©, which is not enough to grow a stel- 
lar mass BH to 10 3 Mq in about 10 s years. On the other 

<C 2 

hand, very light fermions i.e. TO/ ~ 12keV/c would gen- 
erate very massive degenerate FBs of masses greater than 
10 6 Af Q and the BH would grow to 10 9 Af Q in a very short 
time, i.e z « 6.41. 

The growth of seed BHs from DM accretion is investi- 
gated using the quantum cascade mechanism upon which 
low angular momentum DM particles at the inner Fermi 
surface are first consumed by the BH and then due to 
a high degeneracy pressure, higher angular momentum 
particles are pushed inwards and the process continues 
until the entire degenerate FB is consumed by the BH. 
Moreover, molecular clouds have been used as perturbers 
of DM particle orbits outside the FB and we have shown 
that the BH grows faster than the FB. After the BH has 
consumed the entire FB, it then grows by Eddington - 
limited baryonic accretion to higher masses of ~ 10 9 Af© 
at redshifts z ~ 6.41. We also point out that molecular 
clouds of mass 1O 1O M0 have also been detected in the 
host galaxy of the same quasar at a redshift of z ~ 6.41 
(Walter et al. 1500$ . 

We have also constrained the possible starting time of 
accretion, i.e. the time of BH seed formation. From our 
analysis, the mass of a 3 x 10 9 Af Q BH in the quasar SDSS 
J114816. 64+525150.3 at a redshift of z = 6.41 can be 
fitted exactly if the accretion process starts at a time of 
about 2 x 10 8 years, which corresponds to the reionization 
time. 

The seed BH mass is found to be in the range from a 
few solar masses up to an upper limit of ~ 10 4 Af©. For 
a seed BH mass of 10 3 ~ 4 Af Q , Eddington baryonic matter 
accretion would be enough to cause the seed BH to grow 
into a supermassive BH of 3 x 10 9 Af Q mass. The data point 
at a redshift of z = 6.41 can be fitted by only Eddington 
baryonic matter accretion with an efficiency of e ~ 0.01. 

Our model provides a method to find the DM particles 
mass. If it is found that there is a clear lower mass cut 
of 10 3 to lO 6 Af in the distribution of BH masses, then 
this mass can be used for the mass of the FB to find the 
corresponding mass m/ of the fermions which will be in 
the range of 12 keV to 450 keV. If on the other hand the 
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BH mass distribution is a continuous function, then our 
model of BH growth with DM will probably be ruled out. 

The postulated DM particles in this paper were non - 
relativistic at the decoupling time and are usually called 
cold dark matter particles (CDM). The latter have to be 
neutral, stable or quasi-stable and have to weakly inter- 
act with ordinary matter. As mentioned in section 2.1 , 
these particles could be axions which have been investi- 
gated by DAMA/NAI (Bernabei et al. I2()()lj> . The heavier 
particles of mass above lGeV/c 2 could also be of the class 
of DM candidates named WIMPS (Weakly Interactive 
Massive particles). However, in the standard model of par- 
ticle physics, CDM cannot be suitable candidates for par- 
ticles. Thus, a new window beyond the standard model 
of particle physics has to accommodate these particles for 
our model of BH growth to work. The DAMA/NAI ex- 
periment (see Bernabei et al. 2005 for a review) which 
aims at the verification of the presence of DM particles 
in the Galactic halo, will be able to confirm whether GeV 
WIMPS or axions of masses of 12 to 450 keV could exist 
in nature. 

Observations have shown that the masses of supermas- 
sive BHs at galactic centers correlate with the masses of 
the host bulges, i.e. Mbh ~ 0.002 Mbuige (Haering & Rix 
2004). This result is obtained in our model as long as the 
Eddington limited accretion dominates the final growth of 
the BH. This happens for BH masses of more than about 
10 5 M Q (Wang, Biermann & Wandel 12777117)) . 

Mergers of dwarf galaxies as well as the spinning of 
BHs would play an important role in the growth of the 
BHs. The consideration of these two effects will be the 
subject of further investigations. In addition, it would be 
of great interest to study the growth of BHs from boson 
DM particles. Although it has been shown that bosons 
could provide a good fit to the rotation curves in dwarf 
galaxies, it is not yet clear whether an analogous mecha- 
nism could work in galaxies with a 1 /r 2 density fall off. 

While it takes only 8.4 x 10 8 years to grow super- 
massive BHs in most distant quasars, the Galactic center 
might have grown to its current mass of ~ 10 6 M Q with 
only DM accretion in a Hubble time. In the following pa- 
per, we will address the growth of the Galactic center BH. 
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